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Abstract 

The well known Haldane map from spin chains into the 0(3) non linear sigma model 
is generalized to the case of spin ladders. This map allows us to explain the different 
qualitative behaviour between even and odd ladders, exactly in the same way it explains 
the difference between integer and half-integer spin chains. Namely, for even ladders the 
topological term in the sigma model action is absent, while for odd ladders the parameter, 
which multiplies the topological term, is equal to 27rS', where S is the spin of the ladder. 
Hence even ladders should have a dynamically generated spin gap, while odd ladders with 
half-integer spin should stay gapless, and physically equivalent to a perturbed SU (2)i Wess- 
Zumino -Witten model in the infrared regime. We also derive some consequences from 
the dependence of the sigma model coupling constant on the ladder Heisenberg couplings 
constants. 
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1 Introduction 



One of the most studied field theories in 2 dimensions are the non-hnear sigma models. From a 
Particle Physics point of view these models are ideal analogs of 4 dimensional quantum chromo- 
dynamics, since they display asymptotic freedom behaviour [||], dynamical mass generation, and 
existence of instantons 0. In string theory the conformal invariant sigma models are crucial 
to understand the on-shell properties on the strings 0. In Solid State Physics the 0(3) sigma 
model plays also an important role in understanding the properties of spin systems in various 
dimensions. The map from spin chains into the sigma model lead Haldane to the celebrated 
gap of antiferromagnetic Heisenberg chains with integer values of the spin The crucial 

parameter which controls the behaviour of the sigma model is the angle 6 that multiplies a topo- 
logical term into the action and which according to the map takes the value 6 = 2i[S, where 5* 
is the spin of the chain. Haldane's prediction followed from the fact that the sigma models with 
6 = 0(mod27r) are massive field theories This map was deduced in the semi classical 

limit where S >> 1, but there is by now a clear experimental and theoretical evidence of the 
existence of the gap 0,||lO[- For half-integer values of the spin the prediction, based on the 
gap less character of the spin 1/2 chain [0, was that all these models should also be gap less. 
This has been confirmed by numerical computations [0. The sigma model at 6* = vr has also 
been proved to be massless [|T^ . 

Nowadays there is a better understanding of the sigma model at ^ = vr by means of the 
powerful techniques of bosonization |jT^, ||15|, conformal field theory and also from the factorized 
scattering theory []TB[,[0. It has been shown that the low energy physics of the 6 = tt model 
is well described by the SU{2)i Wess Zumino Witten model |jl8[, [|19|. The asymptotically free 
theory in the ultraviolet region, which is described by two goldstone modes, becomes in the 
infrared region a SU{2)i WZW model [0]. This RG flow satisfies the Zamolodchikov c-theorem 
pO| . An alternative description of the sigma model for low energies is given by a marginal- 
irrelevant perturbation of the WZW model by the product of the two chiral SU{2) currents 



18| . In this manner one can do perturbative calculations around the conformal point |21[. It 



is interesting to observe that there are not relevant perturbations of the SU{2)i WZW model, 
which explains why this conformal field theory characterizes the universality class of a large 
variety of spin systems. 

We shall show in this paper that the sigma model methods can be extended to spin ladders, 
which will then allow us to consider certain questions arising in this subject. Spin ladders are 
arrangements of parallel spin chains with nearest neighbour Heisenberg couplings between the 
spins along and across the chains. These spin systems are interpolating structures between 1 and 
2 dimensions. The interest on spin ladders increased enormously when experimentalist discovered 
materials like {VO)2P207 ||2^ and Srn-iCun+i02n [ID, whose magnetic and electronic structure 
was analysized in ^5[. Hence from an experimental and theoretical point of view the spin ladders 
have become a place where to test different ideas concerning strongly correlated systems ( see 
reference I^Gj for a review on the subject). 

A central question in the study of spin ladders is their different qualitative behaviour as a 
function of the number of legs n£. The main conjecture is that ladders with even number of legs 
have a finite spin gap and short distance correlations while the odd ladders have gapless spin 
excitations, and power-law correlations. Many authors have contributed to clarify this question, 
and despite of some initial controversies, it is clear by now that it must be correct |2^-[Q. In 
this paper we shall give further support of this conjecture using sigma model techniques, which 
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will clearly show the topological nature of the mechanism underlying the existence or absence of 
a spin gap in the ladder's spectrum. Our proof is an extension of the Haldane's result from the 
chain to the ladder. Indeed, we shall show that the low lying modes of the ladder are described 
by a sigma model with a value of the angle 9 equal to zero for the even chains and equal to 2ttS 
for the odd ones. Hence from this result and the well known properties of the sigma model at 
= and TT we prove the above conjecture. This kind of proof has been suggested in [34] on 
the basis of the 2d formulation of the sigma model due to Haldane , who showed the absence of 
the topological term in the 2d Heisenberg antiferromagnet [0. However our approach to the 
problem is not really 2 dimensional since we take into account the specific nature of the ladders 
(i.e. objects in between Id and 2d). This different treatment is made clear by the nature of the 
sigma model variables that we use which are one dimensional fields. 

In the case of odd ladders with half-integer spins there is an extension of the Lieb-Mattis- 
Schultz theorem due to Affleck which states that in the infinite length limit , either the ground 
state is degenerate or else there are gapless excitations |^ , . Our results, together with those 
of references [|^,|^,|3y], [^], confirm that the last possibility is the one realized by the 
spin ladders, the spin chain being the particular case = 1. The LMSA theorem works under 
very general circunstances, which is a manifestation of the topological nature of odd ladders and 
half-integer spins. To make this more transparent we shall consider spin ladders with arbitrary 
values of the interchain and intrachain coupling constants. Our results concerning the nature 
of the gap will be independent of the precise values taken by these parameters. We want to 
mention here another topological interpretation of the difference between even and odd ladders 
given in |^ in the framework of the RVB picture According to |^ the even ladders are 



short range RVB systems which have a gap due to the confinement of topological deflects, while 
the odd ladders are long range RVB systems with no confinement and consequently no gap. It 
would be interesting to analize the relation between these two topological interpretations. 



2 AF Spin ladders: goldstone modes 

The Hamiltonian of a spin ladder with ni legs of length N is given by, 

-f^ladder = -f^leg + -f^rung (1) 
^^leg = Eail En=l J a ^a{,n) ■ Sa{n + 1) 
-f^rung = Sa=l ^n=l J'a,a+l Sa('^) ' Sa+i(n) 

where Sa(n) are spin-S matrices located in the a^^ leg at the position n = 1, . . . , A^. We consider 
periodic boundary conditions along the legs ( Sa{n) = Sa{n + N)). The only condition we 
shall impose on the coupling constants Ja and is that they are positive, which guarantee 

that -/^ladder posseses, in the classical limit ( 5* —>■ oo), a minima given by the antiferromagnetic 
vacuum solution, 

Sf^^(n) = (-1)'^+'^ (2) 

where z is the unit vector in the vertical direction. The solution breaks the 0(3) rotational 
invariance of i^iadder down to the subgroup 0(2) of rotations around the z— axis. Consequently 
there should appear two goldstone modes associated to the broken generators and S"^ In 
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the thermodynamic hmit where N oo, with ni keept fixed, the spectrum of the Hamiltonian 
(0) becomes essentially one dimensional, despite of its 2d origin, and hence one expects that 
the quantum corrections will restore the 0(3) symmetry, as it happens for the usual Heisenberg 
chain. Our strategy will then parallel the one used for the study of spin chains. We shall only 
consider the massless degrees of freedom associated to the two goldstone modes and later on we 
shall consider their interaction in the framework of the sigma model. An direct way to find the 
goldstone modes, which are nothing but spin waves, is through the linearized approximation of 



the equation of motion of the spins [40 



The evolution equation of the spin operators of the ladder are given by. 



= i [i/ladder, Sa{n)] = 

-Sa{n) X U {Sa{n + 1) + Sa{n - 1)) + j;,+iS,+iH + J'^^^_,Sa-i{n) 



(3) 



This equation is valid for any a = 1, . . . , with the convenium Jq = J'n^^m+i = ^ind = J^^„. 
Expanding Sa{n) around the classical solution (^), 



one gets in the linearized approximation. 



T class , 



n) + Sa{n) 



(4) 



dCa jn) 
dt 



li-l 



\a+n+l 



s 



Ja {Un + 1) + Un - 1) + 2Ca(ri)) + E^^^ Un) 



(5) 



where Ca(^) = + '^s^^in) and the matrix K^^^ together with a matrix K^^^, which we shall 

use later on, are defined as follows. 



K. 



a,b 



J'a,a+l + J'a,a-l a = ^ 



±J', \a-h\ = l 



(6) 



The solution of eqs.(||) are given by plane waves, 
Introducting (|^) into one gets. 



(7) 



^Mk) = 4sin^ f Ja Mk) + EbKb Mk) 



^Uk) = 4cos2 f J„ ^aik) + EbKb Mk) 

These equations have massless and massive modes in the limit — * oo ( periodicity along 
the legs implies k = 27rm/N ,m = 0,1, ■ ■ ■ , N — 1). These modes can be obtained by expanding 
uj{k), ipa{k) and (paik) in powers of the momenta k. For the massless modes this expansion reads. 



u= vk + 0{k^) 
Mk)= kAa + 0{k') 

Mk)= Ba + eCa + 0{k^) 



(9) 
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The equations for the coefficients Aa, Ba and Ca follows from (j^), 



0= nK-^Bb (10) 

^Ba= Y.bLa,bA, (11) 

JaBa+TbK^bCb (12) 

where we have introduced yet another matrix La^b given by, 

La,b = ^Ja Sa,b + (13) 

which is going to play an important role in the construction. Notice that La,b is a positive 
definite matrix. 



The solution of equation (IC) is uniquely given hj Ba = B Va. This result follows from the 
observation that K~f^ is a generalized incidence matrix associated to a graph consisting of 
points labelled by a and links joining the points a and b whenever ^ is non null. Since the 
graph is connected (i.e. ^ is a non singular matrix) there is a unique vector satisfying equation 
(p!0|). Connectedness of the graph simply means that the ladder cannot be split into two or more 
subladders. This graph together with its incidence matrix contains all the information of the 
rungs of the spin ladder relevant to our problem. 



Aa = ^BY.L-i (14) 
b 



The solution of (|TT]) is given by, 



where we have inverted the matrix L ( recall that L is positive definite). To solve equation (|12D 
we proceed in two steps. First of all we sum over the index a in (|T2|) and use the fact that 
J^a^ab — 0; to of term proportional to Ca- This give us an equation for the spin 

wave velocity v, which with the help of (0) can be written as, 

VV T.aJa ^^^^ 



-S ' Ha,bLa^b 

Both the numerator and the denominator of (^) are positive which yields a real value for 
v/S. 



The solution of ( |12D for the vector Ca constitute in fact a one parameter family of solutions 
given by Ca + xBa with x arbitrary. This freedom refiects the linearity of equations (P). Multi- 
plying the whole solution (^ by a A; dependent factor produces a term of the form k^Ba in (^. 
The "transverse" components of Ca can then be obtained by inverting the matrix K~ matrix in 
the subspace orthogonal to its zero eigenvector. 

Next we shall briefiy consider the massive modes. The value of the gap A = ci;(/c = 0) can 
be simply obtained setting = in (R|), 



^M^) = T.bK-M^) (16) 

f0a(O) =E6i^a,6V^6(0) 

Hence combining both equations we get that (A/S")^ is given by the eigenvalues of the matrix 
LK~ or alternatively K~L. One of this eigenvalues is zero and corresponds to the massless mode 
studied above and the others are all non zero and positive corresponding to the massive modes. 
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3 cr— model mapping 



Let us recall how one maps the Heisenberg spin chain into the Id cr— model (we shall follow 
closely reference P). The spin wave analysis shows that the spin operators S(n) has two smooth 
components centered and momenta k = and k = tt, which can be identified with the total 
angular momenta 1 and the staggered field ip respectively. The relation between these operators 
can be written as follows, 



S(2n) = - Sip{x) (17) 
S(2n + 1)= l{x) + S^{x) 

where x = 2?7, + | is the midpoint coordinate of the block formed by the points 2n and 2n + 1. 
Inverting eq.(|T7D one gets, 



l{x)= (S(2n + l) + S(2n))/2 
(^(x) = (S(2n + 1) - S(2n))/25 

1 and if satisfy the following commutation relations. 



:i8) 



r(x),P(y) = 
r{x),^^{y) = 



(19) 



which can be derived from the SU{2) algebra satisfied by the spin operators S{n). The term 
is the lattice version of the Dirac's delta function 5{x — y). The factor two in the denominator 
is simply the lattice spacing of the 2-block arrangement of the chain. The fact that S(n) are 
spin-S matrices satisfying the relation S^(n) = S{S + 1), imply two additional equations for I 
and (p, namely. 



= 1-17^2 + 0(1/^) ^ 1 



(20) 



Equations (0) and (pOD are the standard relations satisfied by the a— field (f and the angular 
momenta 1. Introducing eqs.(p!7D into the spin chain Heisenberg Hamiltonian one gets, taking the 
continuum limit and eliminating higher derivatives terms, the standard a— model Hamiltonian, 



2 



dx 



9 



(21) 



where y?' = d^f. The theta parameter, coupling constant and spin velocity take the following 
values, 



e = 2nS, g=-, v^ = 2JS (n, = 1). 



(22) 



5 



The Hamiltonian (21) can be obtained from the 2d cr— model Lagrangian, 



L = + — e"*' V ■ (d^V X d,^) (23) 



1 

It is our aim to generahze the previous construction to spin ladders. First of all we shall 
divide the ladder into blocks of two consecutive rungs and define a— model variables for each of 
them. The spin wave analysis of the previous section suggest the following ansatz, 



Sa{2n)= Aa\ + la + S{-iy{y^ + y^a) 

S,(2n + 1)= AJ + la + Si-ir+'iif + ifa) 



(24) 



where I and (f are the candidates for the a variables and la and (fa are some extra slowly varying 



fields needed to match the number of degrees of freedom in both sides of (24). For these to be 
the case we shall impose the following "transversality conditions" on 1^ and (fa, 



J2a Va 





: 



(25) 



which are only needed for nf > 1. 
operators as follows, 



Using 



we can express 1 and cp in terms of the spin 



K^) = Ea [Sa(2n + 1) + Sa{2n)] / (2 A,) 
^(x) = Eai-ir^' [Sa(2n + 1) - Sa(2n)] / (25n,) 



(26) 



We want 1 and ip to satisfy the algebraic relations (19), which can be achieved imposing. 



A. 



Lai 



where we have used eq.(|4D. Similarly I and f as given by (p6D satisfy eqs. similar to 



(27) 



= 1 + OiXjSne) (28) 
I-<^ = 0(l/5n,) 

Hence in the limit Sn^ >> 1 we obtain the constraints which define the a— model. From 
(^) it seems that the expansion parameter that we are employing is Sni rather than S. If this 
is correct then considering higher spins is equivalent, from the sigma model point of view, to 
considering ladders with many legs. We shall return later on to this suggestion. Another point 
which is worth to mention is that \[x) and fix) represent total angular momenta and staggered 
magnetization of the rung taken as a whole. This is why they are Id densities depending only on 
the single coordinate x. Given the relations ( P^ we can now write the spin ladder hamiltonian 
(0) in the variables 1, and as follows. 
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Hladder — I ^ {j2a,b La,b (^a^fe 1^ + L h) 

Ea Ja + V'af + Ea,6 K~, ip^ ifb (29) 
+2SY.a{-lYJa [(Aal + L) + ^'J + + V'a) (^al + L)] } 

To derive (^Up we have used ([TT|), (P3| ) and the following formula, 

(AJ + \af + S''{ip + ifaf = S{S + 1) (30) 
(A J + 1,) {if + ^a)=^ 

which is a consequence of the relations S^(n) = S'(S'+1) (since we are working in the semiclassical 
limit S* >> 1 we shall keep only the highest power in S). To decouple the fields (p and (pa in (p9D 
we have to choose the same value of Ja for all the legs. Indeed upon this condition the cross term 
Ea Ja^'^'a (p9|) vauishes as a consequence of (pSj). We thus obtain that 99 is a massless field 
while the fields are massive. Let us now concentrate on the massless field, whose Hamiltonian 
reads, 



(massless) _ f dx 
ladder 



/ f [(Ea,fe La,,AaA,) F + JaV'' + 25 Ea(-l) (1 + 1)] (31) 



This is precisely the a— model Hamiltonian given in ( PT| ) with an appropiate identification 
of 6,g and v. Let us first consider 6 whose values is given by. 



A simplification of (|3^) is achieved using (p^), 

= 87r5 5:(-l)'^+V„L-^ 



(32) 



(33) 



a, 6 



A convenient way of writing ( P^D is by means of two n^-dimensional vectors \F > and \AF > 
defined as follows, 



\F >-- 



whose scalar product is. 



/ 1 \ 



V 1 / 



< AF\F >-- 



\AF >-- 



( 1 \ 



(34) 



rii : even 



Using 



Til : odd 

we write (031) in matrix notation as, 



(35) 
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e = SnSne < AF\ J L"^ \F > (36) 
where J is a diagonal matrix whose entries are Ja- Recalhng the well known operator identity, 



B 



we transform (|3q) into, 



A + B A A A + B 



9 = 2nSne < AFl 1 - K 



\F > 



(37) 



(38) 



4J + K+, 

Then noticing that the vector | AF > is annihilated by and using eq. (p5D we arrive finally 



at , 



111 : even 



9 



(39) 



2713 ne : odd 

Taking into account that 9 is defined modulo 2n then we can write eq.(p9D simply as, 

9 = 2nSni (40) 

This result is valid for any value of the coupling constants Ja and Ja,a+i as long as they 
are non vanishing. This confirms the topological nature of the result in total agreement with 
the LMSA theorem. Moreover the derivation of ( ^OD suggest that the 9 term for spin ladders is 
related to the transition amplitude from ferromagnetic to antiferromagnetic configurations along 
the rungs (^5]) . A path integral derivation of (^0]) , along the lines of |^ , would probably throw 
some light on this interpretation. 

Next we shall give the expressions of the a— model coupling constant g and the velocity v^, 
in terms of the ladder parameters. 



S 



a,b,c 



1/2 



2^; — - K 



1 



(2E 



b,c J^b,c 



(41) 



(42) 



where Sn^ is equal to 1 (or 0) whenever ne is odd (or even). 

Compairing eqs.(^) and (0) for = 1 we get that = v, but for ne > 1 the two velocities, 
fg- and V, do not coincide (for ne even one has = \/2 v). We interpret this fact as a kind of 
interference effect between the legs of the ladder which makes the effective spin velocity Uo- to 
differ from the spin wave velocity v. The most interesting case for practical applications is when 



Ja = J and J' 



a,a+l 



J' Va. We shall give below the values of g and Va in this situation. 



S 



^f{ne, J' /J) - \ 5^, 
4J 

ne g fine, J' / J) 



1/2 



(43) 



(44) 
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The function /(n^, J' / J) appearing in these formulae is defined as, 



/K,JVJ)=<F| (l + ff)"' \F> 



(45) 



^ni, + 2 I]m=l,3,---,n^-l (^^^ 1^ 



1 + 4 C0S2 ^ 



In the particular cases where rii = 2 and 3 we obtain, 



for n£ = 2 



SV2 



1 + 



1/2 



2^257 (l + 



1/2 



(46) 



for n£ 



2 

s 



4,7 



' ^ 4,/ 



1/2 



2JS 
3 



1/2 



(47) 



1+1 



We have assumed in ( ^T] ) that S' is half integer. From eq.(^3|) we derive that g{ne,J'/J) is 
a monotonically increasing function of the ratio J'/ J, which implies that the ladder is more 
disordered in the strong coupling regime that it is for weak coupling. In fact we get. 



lim q 

J'/J^oo 



2 
S 

V2 

S 



{J'/J) 



1/2 



(4^ 



oo ni : even 

: odd, 5* : half integer 
ne : odd, S : integer 

This equation shows that the difference between even and odd ladders appears not only in 
the topological term but also in the behaviour of 5^ as a function of the ratio J' / J in the strong 
coupling limit. For ne even the sigma model enters into the strong coupling regime g » 1, 



which is dominated by the angular momentum term P in the sigma model Hamiltonian (pT]). 
Let us suppose that we discretize the sigma model Hamiltonian at 6' = as in references 0, 



2 



E 



g l^(^) '^{n) ■ {p{n + 1) + cte 



(49) 



with \{n) satisfying the standard angular momenta algebra, such that F has the spectrum 
/(/ + 1), Z = 0, 1, . . .. l(n) gives the angular momenta of the n^^ rung of the ladder. Then in the 
strong coupling limit the ground state of (|49|) is obtained choosing the representation / = for 
each n. The first excited state has / = 1 at one site and energy gva, which is the value of the 
gap in the limit g » 1. In the case ni = 2 we get from eqs (^61), 

v^g ~ J', for J'/J » 1 (50) 

The second term in ( ^91) produces shifts in the ground state energy and also delocalizes 
the / = 1 excitation producing a band of states. The gap only vanishes at g = 0. These 



results agree with the ones obtained using very different techniques namely, numerical |27|, |28 
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renormalization group |29| , mean field [3^, finite size and bosonization |35|. However in 



order to claim full agreement we have also to analize what happens with the other massive modes 
that we discarded in the mapping of the ladder Hamiltonian into the sigma model Hamiltonian. 
If the mass which is generated dynamically by non perturbative effects for the field ip is smaller 
than the gap associated to the massive modes (pa then expect that the map must be essentially 
correct, except for a finite renormalization of the coupling constant g and the spin velocity v„ 
This issues will be considered elsewhere. 

For the odd ladders the asymptotic value of g is in agreement with (^2]), in the sense that 
the odd ladders with spin 1/2 can be though of as single chains with a spin 1/2 and an effective 
couphng constant Jeff. Indeed, for S* = 1/2 we get from (0) that g = 4, which is the same value 
we obtain in ( p2|) for the single spin 1/2 chain. 

Let us consider now the weak coupling limit J' / J << 1. From ( ^31) and (|45|) we get. 



v^/ {Srii) Hi : even or ni : odd, S : integer 

lim g = { (51) 

2/ [S^2nj - l) Ui : odd, S : half integer 

Which implies that g depends essentially on the combination Sn£, as we anticipated in the 



discussion of eq. ([281) . The isotropic case , J = J', is in fact closer to the weak coupling values 
(^) than to the strong coupling ones (^Sf) . Thus for Srii >> 1 the value of g will be small 
and we may use the formula exp{—27i/g) to estimate the value of the energy or mass scales of 
the system. This implies in particular that the mass gap for the even ladders with ne large will 
decrease as exp{— cte ne). This agrees at least qualitatively, with the numerical results which 
give a spin gap Agpin at the isotropic case equal to 0.504J for ne = 2 and Agpin ~ 0.2 J for ne = 4. 
Thus in the limit — oo the gap of the even ladders should vanish exponentially. As the odd 
ladders are already gapless for any number of legs one reaches in the limit the same result 
for both even and odd chains. However one must be careful in this limit since as we mentioned 
above the massive modes that we discarded in our mapping to a Id sigma model are becoming 
more important as we increase the number of legs. A more careful analysis of this questions is 
needed. 



4 Final Considerations 

The application of the sigma model techniques to spin ladders has allowed us not only to confirm 
the topological origin of the qualitatively different behaviour of even and odd ladders, but also 
to get some hints about the dependence of the physical quantities on the values of the coupling 
constants J and J'. Much work remains to be done in this direction, but we believe that the 
sigma model offers an unified, economic, and consistent approach to spin ladders. The connection 
we have stablished in this paper allows in principle to apply the knowledge accumulated in the 
past in the study of the sigma model to the understanding of spin ladders. 

An interesting "recent" result concerning the sigma models at = vr is the proof of its exact 
integrability [|I^,|T^ , which is the parallel of the well known integrability of the sigma model 



at ^ = 0. The proof of integrability is done in the framework of the factorized scattering 
theory. For 6 = the S-matrix is formulated for a 0(3)-triplet of massive particles, while for 
6 = n there are two 0(2)-doublets of left and right moving massless particles, whith three 
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types of scatterings: left-left, right-right and left-right. Using the powerful techniques of the 
thermodynamic Bethe ansatz one can compute finite size effects of various observables. In this 
way one can prove that the RG-flow for 6 = n, goes from the UV asymptotically free model 
with c=2 to the IR massless 5'[/(2)iWZW model with c=l, as we indicated in the introduction. 
The results one gets using these exact techniques agree at the perturbative level with the ones 
obtained from the perturbation of the WZW model by the marginal irrelevant operator JlJr 
|19| |. In this way one explains the logarithmic departure from the scale invariant results which 
can then be compared with experimental or numerical results . These logarithmic corrections 
all depend on a mass scale A, which is generated dynamically in the sigma model, and which 
for small values of g is given essentially by 1/a exp{—2Ti/g), with a the lattice spacing. In the 
factorized S-matrix theory the parameter A appears explicitly in the expression of the energy 
and momentum of the particles. An important problem is to derive the relation between A and 



the microscopic parameters of the model appearing in the Hamiltonian, namely n^, J and J' |^ 



If the Hamiltonian happens to be integrable then one should be able to find an exact expression 
for A, but in general this will not be possible and so one has to use some approximation method. 
Numerical computations of thermodynamics quantities of the spin ladders and they comparison 
with the theoretical predictions may also be very useful in stablishing this connection 
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